Effect of electron and hole doping on the structure of C, Si, and S nanowires 
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We use ab initio density functional calculations to study the effect of electron and hole doping on 
the equilibrium geometry and electronic structure of C, Si, and S monatomic wires. Independent 
of doping, all these nanowires are found to be metallic. In absence of doping, C wires are straight, 
whereas Si and S wires display a zigzag structure. Besides two preferred bond angles of 60° and 
120° in Si wires, we find an additional metastable bond angle of 90° in S wires. The equilibrium 
geometry and electronic structure of these nanowires is shown to change drastically upon electron 
and hole doping. 

PACS numbers: 73.22.-f, 61.46.-w 



The close relationship between geometric and elec- 
tronic structure of molecules and solids underlies the rich- 
ness of phenomena described by Physics and Chemistry. 
The electronic structure plays a central role in determin- 
ing the equilibrium geometry by minimizing the corre- 
sponding free energy. Geometrical structure, in turn, de- 
termines the electrostatic field confining the electrons. 
This subtle two-way relationship between geometry and 
electronic structure underlies the complex behavior of 
even simple systems, such as monatomic nanowires or 
triatomic molecules. In the latter case, rather heuris- 
tic rules, such as the Walsh rulei, have been used to 
rationalize, whether a trimer carrying a particular num- 
ber of electrons should be linear or bent. Approaches 
based on the theory of directed valence, including the 
Valence Shell Electron Pair Repulsion (VSEPR) model^, 
have generally proven useful in determining the bonding 
geometry in polyatomic molecules. These approaches, 
however, miss subtle differences in bonding between el- 
ements of the same group and are too rough to predict 
the effect of doping on the bonding geometry. 

Most previous studies have focussed on elements with 
a ductile, metallic bulk phase due to their capability 
to form nanowires by mechanical stretching. Interest 
in the geometry of such nanowires has been triggered 
by electron microscopy observations of monatomic Au 
nanowires with atomic resolution, and the correlation be- 
tween their atomic structure and quantum conductance^ 
Even the structure of the simplest nanowires is still a 
matter of debate, since the observed inter-atomic dis- 
tance was found to be significantly larger than the bulk 
inter-atomic distances. Several explanations have been 
offered as an interpretation of the initial experimental 
data. These include a rotating zigzag wire structure^ 
stabilization of the nanowire by impurity atoms with 
a low cross-section for electrons*^ and charging effects 
that would expand the inter-atomic distance. 7 

Inspired by these investigations on Au, theoretical 
studies have been expanded to nanowires of other met- 
als, including Al, Ag, Pd, Rh, and Ru. 8 Some of these 
nanowires are reported to exhibit zigzag deformations 
with more than one stable bond angle at zero strain. In 
contrast to VSEPR model predictions, theoretical stud- 



ies of group IV elements suggest that C forms a linear 
wire, whereas the Si forms a zigzag wire^i^ Nanowires 
of non-metals are of interest as either free-standing sys- 
tems, or constituents of molecular electronics devices^ 
In the latter case, their geometrical and electronic struc- 
ture, and thus their transport properties, may be influ- 
enced by charge doping caused by the chemical environ- 
ment and the presence of a current. 

In the present study, we use ab initio density functional 
calculations to investigate the effect of charge doping and 
structural constraints on the equilibrium geometry and 
electronic structure of C, Si, and S monatomic wires. For 
different levels of doping, we present total energy plots 
as a function of lattice constant and bond angle, which 
provide important information about the structure and 
stiffness of these nanowires, when exposed to particular 
environments or deformations. We find neutral C wires 
to be straight, whereas Si and S wires display a zigzag 
structure. Besides two preferred bond angles of 60° and 
120° in Si wires, we find an additional metastable bond 
angle of 90° in S wires. We also observe changes in the 
equilibrium geometry and electronic structure induced by 
doping, including a zigzag distortion of electron doped 
C nanowires and the number of stable geometries in S 
nanowires depending on the doping level. 

Our calculations for C, Si, and S monatomic chains 
are based on the density functional theory (DFT) 12 : 13 
within the local density approximation (LDA). We use 
Troullier-Martins ab initio pseudopotentials to describe 
the interaction of valence electrons with atomic cores^ 
and the Perdew-Zunger form of the exchange-correlation 
potentiali^^, as implemented in the SIESTA code-i 7 -^ 
Our basis consists of doublc-zeta localized orbitals with 
polarization functions (DZP). The range of the local- 
ized orbitals is limited in such a way that the energy 
shift caused by their spatial confinement is no more than 
100 meV i 19 ' 20 i 21 The charge density and potentials are 
calculated on a real-space grid with a mesh cutoff energy 
of 200 Ry for C and 150 Ry for Si and S. This is sufficient 
to achieve a total energy convergence of 1 meV/atom. 
We use periodic boundary conditions in all directions 
and separate the nanowires laterally by 9 A to prevent 
them from interacting. We sample the ID Brillouin zone 
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FIG. 1: (Color online) Schematic of a unit cell of (a) a single- 
stranded and (b) double-stranded nanowire, with the unit cell 
of length a delimited by the dashed lines. The key distinguish- 
ing feature between the two structures is the value of the bond 
angle 8. The atomic positions are indicated by solid circles, 
nearest neighbor bonds by solid lines, second neighbor bonds 
by dashed lines. 



of the nanowires by 96 k-points for all the lattice con- 
stants considered. Since electronic structure calculations 
in a superlattice geometry require global charge neutral- 
ity, we perform calculations for charged chains on the 
background of a uniformly distributed counter-charge. 

For a given lattice constant, we globally optimize each 
system to determine its optimum geometry, total energy, 
and electronic structure. We accommodate two atoms 
per unit cell, and use an initial geometry with two equiv- 
alent bonds as a starting point of the optimization. We 
identify each nanowire structure by the pair of values 
(a, 8), where a is the lattice constant in Angstroms and 9 
the bond angle in degrees, as depicted in Fig.QJa). Obvi- 
ously, for bond angle values 9<60°, the second neighbor- 
interaction along the wire direction may dominate over 
the nearest neighbor interaction, causing a spontaneous 
transition to a double-chain. This geometry is depicted 
in Fig. QIb). We consider a structure to be optimized 
when none of the residual forces acting on atoms exceed 
the value of 0.01 eV/A. 

Fig.[5]shows contour maps of the total energy of carbon 
nanowires £' tot (a, 6) as a function of the lattice constant 
a and bond angle 6. Independent of doping, we generally 
find the bond angle to increase with increasing lattice 
constant. Especially in systems where the bond length 
is stretched beyond the equilibrium value, this trend re- 
flects the tendency to optimize the bond length first and 
only then the bond angle, suggesting that bond stretch- 
ing is harder than bond bending. 

For a neutral carbon wire, depicted in Fig. [DJa), the 
equilibrium structure is found at (2.6,180), correspond- 
ing to a linear chain. When the bond angle is con- 
strained to 0<6O°, a secondary minimum evolves at 
(1.4,20), corresponding to a double-chain with only a 
slightly longer nearest neighbor distance in each strand. 
For a hole-doped wire carrying +0.5 electrons/atom, de- 
picted in Fig. [2b), the optimum structure is still lin- 
ear and the interatomic spacing is unchanged, corre- 
sponding to (2.6,180). In contract to the undoped sys- 
tem, however, the Coulomb repulsion suppresses the sec- 
ondary minimum corresponding to a paired wire config- 
uration. Upon electron doping at the level of —0.5 elec- 
trons/atom, depicted in Fig. (He), the global geometry 
optimum shows a much weaker dependence on the bond 
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FIG. 2: (Color online) Energy of a C nanowire as a function 
of the lattice constant a and the bond angle, for an undoped 
system (a), and for doping levels of +0.5 electrons/atom (b), 
and —0.5 electrons/atom (c). Contour lines, representing the 
total energy per unit cell, are separated by 0.2 eV. Optimized 
geometries for each lattice constant are shown by the white 
dotted lines. Lines of constant bond length are indicated by 
black solid lines. The electronic density of states for the glob- 
ally optimized geometries, denoted by the solid x in (a-c), are 
shown in (d-f), respectively. The Fermi level lies at E — 0. 



angle than in neutral or hole doped wires. The sys- 
tem gains ~30 meV/atom, when relaxing from the linear 
structure at (2.7, 180) to a bent geometry, characterized 
by (2.6, 150). Increased electron doping up to —0.7 elec- 
trons/atom further changes the bond angle to 6*«140° 
and the bond length to 1.47 A at a still larger energy 
gain with respect to the linear structure with 1.3 A as 
optimum bond length. Also in electron doped carbon 
nanowires, the secondary minimum corresponding to a 
paired wire configuration is suppressed. Even though 
the interaction between neutral wires is negligible at the 
inter-wire separation of 9 A, the Coulomb repulsion be- 
tween charged wires is strong enough to suppress the sec- 
ondary minima associate with paired wires. This has 
been confirmed by observing a stabilization of charged 
nanowire systems when increasing the inter-wire separa- 
tion. 

The robustness of the linear structure of carbon chains 
is also reflected in the electronic density of states for the 
globally optimized geometries. As seen in Figs. [2f d-f), 
the density of states is rather featureless, suggesting that 
electron or hole doping does not modify the nature of 
electronic states responsible for bonding. Whether neu- 
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FIG. 3: (Color online) Energy of a Si nanowire as a function 
of the lattice constant a and the bond angle, for an undoped 
system (a), and for doping levels of +0.5 electrons/atom (b), 
and —0.5 electrons/atom (c). Contour lines, representing the 
total energy per unit cell, are separated by 0.2 eV. Optimized 
geometries for each lattice constant are shown by the white 
dotted lines. Lines of constant bond length are indicated by 
black solid lines. The electronic density of states for the glob- 
ally optimized geometries, denoted by the solid x in (a-c), are 
shown in (d-f), respectively. The Fermi level lies at E — 0. 



FIG. 4: (Color online) Energy of an S nanowire as a function 
of the lattice constant a and the bond angle, for an undoped 
system (a), and for doping levels of +0.5 electrons/atom (b), 
and —0.5 electrons/atom (c). Contour lines, representing the 
total energy per unit cell, are separated by 0.2 eV. Optimized 
geometries for each lattice constant are shown by the white 
dotted lines. Lines of constant bond length are indicated by 
black solid lines. The electronic density of states for the glob- 
ally optimized geometries, denoted by the solid x in (a-c), are 
shown in (d-f), respectively. The Fermi level lies at E — 0. 



tral or doped, carbon nanowires arc metallic, with a 
nearly constant density of states near the Fermi level. 
The van Hove singularities in the spectrum are charac- 
teristic of one-dimensional systems. 

Contour maps of the total energy E to t{a,0) of silicon 
nanowires are shown in Fig. [3l The potential energy 
surface of undoped Si nanowires exhibits two minima at 
(3.8, 120) and (2.4,60), as shown in Fig.[^a). The bond 
length tf«2.2 A in the structure with 6*~120° is shorter 
than in the double-stranded structure, formed by equilat- 
eral triangles with d«2.5 A. As seen in Figs. |3jb-c), the 
shape of the energy surface changes only slightly upon 
doping. In particular, the two stable geometries remain- 
ing nearly unchanged, but the one with the smaller bond 
angle becomes more stable upon doping. 

Even though both silicon and carbon are group IV el- 
ements, the geometries of these nanowires are very dif- 
ferent. The underlying cause is the different electronic 
structure of these nanowires, reflected in the electronic 
density of states for the globally optimized geometries. 
Same as C chains, also Si chains maintain a metallic char- 
acter. In contrast to Figs. [2j d-f), the density of states of 
Si nanowires in Figs. [3Jd-f) is very rich in features, sug- 



gesting that the equilibrium geometry is largely deter- 
mined by the energetic availability of particular orbitals 
for bonding. 

As suggested previously^ and seen in Fig.fJJa), the po- 
tential energy surface of S nanowires is more structured 
than that of C and Si nanowires. The global minimum in 
the neutral system is at a = 3.6 A and 9 = 120°. Similar 
to C and Si nanowires, we find a second shallow minimum 
near 9 = 60°. In addition to the 9 = 120° and 60° bond 
angles, found in other nanowires, the S wire shows an 
additional minimum at 9 — 90°. Not shown in Fig. 0Ja) 
is the fourth minimum at (2.2,40), corresponding to a 
pair of linear wires with the bond length d = 2.2 A. The 
bond length in the zigzag wires with 9 — 90° and 120° 
is d = 2.1 A and increases to 2.4 A at 6 = 60°. In con- 
trast to Si wires, doping by electrons and holes changes 
the total energy surface of S wires drastically. The sys- 
tem carrying a net charge +0.5 electrons/atom, depicted 
in Fig. Bfb), has only one minimum at a = 3.0 A and 
9 = 90°. There is also an indication of a second spurious 
minimum at a = 2.2 A and 9 < 60°, indicating the pref- 
erence of the system in a short unit cell to dissociate into 
two parallel wires that repel each other. In the electron- 
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FIG. 5: Total energy changes in neutral and doped monatomic 
nanowires of C, Si and S as a function of the imposed lattice 
constant a, with the bond angle 6(a) given by the white dotted 
lines in Fig.[2Ja-c) for C, Fig-EHa-c) for Si, and Fig.@|a-c) for 
S. 



doped system at the level of —0.5 electrons/atom, de- 
picted in Fig. BJc), we observe two minima at (4.2, 135) 
and (3.0,80). There is also a spurious additional mini- 
mum near (2.2,27), an indication of the instability of a 
charged S nanowire, crowded in a short unit cell, with 
respect to two parallel wires subject to Coulomb repul- 
sion. 

The electronic density of states of S chains, depicted 
in Figs. H)(d-f), is very rich in features, same as that 
of Si nanowires in Figs. [3jd-f ) , thus providing a back- 
ground for explaining the number of stable morphologies 
and their deviation from heuristic models. Similar to C 
and Si chains, the electronic spectrum of sulfur nanowires 
shows no gaps near the Fermi level, indicating metallic 
behavior. 

Total energy differences in neutral and doped C, Si and 
S nanowires are shown in Fig. [5] as a function of the lat- 
tice constant a along the optimum trajectories 0(a), given 
by the white dotted lines in Figs. H]|3Ja-c). Clearly visi- 
ble is the number of energy minima, which changes from 
element to element. A second, clear message emerges, 
namely that electron or hole doping plays a much less 
significant role in carbon nanowires than in those of sili- 
con and, even more so, of sulfur. 

For the charge neutral systems, the equilibrium bond 
angles in nanowires follow the trends found in trimers. 



Both the C trimer and nanowire are linear, whereas the 
Si trimer and nanowire are bent, with the bond angles 
of ss80° in the trimer and the secondary minimum at 
«120° in the Si nanowire. Similar to the S nanowire, the 
S trimer exhibits a more complex energy surface with two 
minima at 9 = 120° and 60°. The structural difference 
between C3 and Si3 had been formerly partly attributed 
to the role of d orbitals in Si3^ Our results differ from 
these findings, since Si3 bends spontaneously even when 
the d orbitals are absent^ 

Even though the initial nanowire geometry, character- 
ized by (a, 6), had two equivalent bonds per unit cell, the 
optimizations have been performed with no further con- 
straints. We found that also all our final structures have 
two equivalent bonds, but cannot exclude the possibility 
of more stable structures with inequivalent bond lengths. 
In the present work, we confined ourselves to finding out, 
whether particular wires are linear or bent. Even though 
bond length alternation caused by a Peierls instability is 
an interesting problem, its exhaustive study exceeds the 
scope of the present work. Reliable results are not easy 
to come by, since calculated Peierls distortions appear to 
partly depend also on the computational approach, with 
distortions predicted by LDA/GGA generally smaller 
than those based on Hartree-Fock calculations i 10 i 25 

In summary, we used ab initio density functional cal- 
culations to investigate the effect of charge doping and 
structural constraints on the equilibrium geometry and 
electronic structure of C, Si, and S monatomic wires. Un- 
like bulk C, Si, and S, the nanowires of these elements are 
found to be metallic. For different levels of electron and 
hole doping, we determined the total energy as a function 
of lattice constant and bond angle. These results provide 
important information about the equilibrium structure 
and stiffness of the nanowires, when exposed to particular 
environments or deformations. We found neutral C wires 
to be straight, whereas Si and S wires displayed a zigzag 
structure. Besides two preferred bond angles of 60° and 
120° in Si wires, we found an additional metastable bond 
angle of 90° in S wires. We found carbon nanowires to 
be linear and rigid, with a tendency to zigzag distortion 
at moderate electron doping. Whereas the zigzag geome- 
try of silicon nanowires is not affected by doping, that of 
sulfur nanowires changes strongly, with even the number 
of stable geometries depending on the doping level. 
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